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Abstract. We introduce here the two-fractal model of earthquake dynamics. As the fractured 
surfaces have self-affine properties, we consider the solid-solid interface of the earth's crust and 
the tectonic plate below as fractal surfaces. The overlap or contact area between the two surfaces 
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, give a measure of the stored elastic energy released during a slip. The overlap between two fractals 
change with time as one moves over the other and we show that the time average of the overlap 
distribution follows a Gutenberg-Richter like power-law, with similar exponent value. 
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Introduction 



i , 

£H , The earth's solid outer crust, about 20 kilometers in average thickness, rests on 
O ' the tectonic shells. Due to the high temperature-pressure phase changes and the 
, consequent powerful convective flow in the earth's mantle (a fluid of very high 
density), at several hundreds of kilometers of depth, the tectonic shell, divided 
^ ■ into a small number (about ten) of mobile plates, has relative velocities of the 
. order of a few centimeters per year (Gutenberg and Richter, 1954; Kostrov and 
£T) ' Das, 1989; Scholz, 1990). Over several tens of years, enormous elastic strains 
develop on the earth's crust when sticking (due to the solid-solid friction) to the 
moving tectonic plate. When sudden slips occur between the crust and the tectonic 
' plate, these stored elastic energies are released in 'bursts', causing the damages 
O . during the earthquakes. 

Earthquakes occur due to fault dynamics in the lithosphere. A geological fault 
is created by a fracture in the rock layers, and is comprised of the rock surfaces in 
contact. The two parts of the fault are in very slow relative motion which causes 



^ ' the surfaces to slide. Because of the uniform motion of the tectonic plates, the 



elastic strain energy stored in a portion of the crust (block), moving with the plate 
q relative to a 'stationary' neighboring part of the crust, can vary only due to the 
random strength of the solid-solid friction between the crust and the plate. A slip 
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occurs when the accumulated stress exceeds the resistance due to the frictional 
force. The potential energy of the strain is thereby released, causing an earthquake. 
As mentioned before, the observed distribution of the elastic energy release in 
various earthquakes seems to follow a power law. 

The slip is eventually stopped by friction and stress starts developing again. 
Strain continues to develop till the fault surfaces again slip. This intermittent stick- 
slip process is the essential characteristic feature of fault dynamics. The overall 
distribution of earthquakes, including main shocks, foreshocks and aftershocks, 
is given by the Gutenberg-Richter law (Gutenberg and Richter, 1954; Gutenberg 
and Richter, 1944): 

log 10 Nr(M >M) = a-bM (1) 

where, Nr(Al > M) denotes the number (or, the frequency) of earthquakes of 
magnitudes M that are greater than a certain value M. The constant a represents 
the total number of earthquakes of all magnitudes: a = log 10 Nr(Al > 0) and the 
value of the coefficient b is presumed to be universal. In an alternative form, the 
Gutenberg-Richter law is expressed as a relation for the number (or, the frequency) 
of earthquakes in which the energy released £ is greater than a certain value E: 

Nr(£ > E) ~ E~ hlli , suggesting Nr(£) ~ E~\ (2) 

for the number density of earthquakes, where y - 1 + b/p. The value of the 
exponent y is generally observed to be around unity (Knopoff, 2000); see also 
http : //web . cz3 . nus . edu . sg/ " chenk/gem2 503_3/notes7_l . htm. 

One class of models for simulating earthquakes is based on the collective 
motion of an assembly of connected elements that are driven slowly, of which 
the block-spring model due to Burridge and Knopoff (Burridge and Knopoff, 
1967) is the prototype. The Burridge-Knopoff model and its variants (Carlson 
et al, 1994; Olami et al., 1992) have the stick-slip dynamics necessary to produce 
earthquakes. The underlying principle in this class of models is self-organized 
criticality (Bak, 1997). 

Another class of models for simulating earthquakes is based on overlapping 
fractals, which will be discusses in details in the next sections. 

2. Fractals 

A fractal is a geometrical object that displays self-similarity on all scales. For 
random fractals, the object need not exhibit exactly the same structure at all scales, 
but the same 'type' of structures appear on all scales. For example, the Black sea 
coastline measured with different length rulers will show differences: the shorter 
the ruler, the longer the 'length' measured but not exactly following the ratio of 
the inverse of the lengths of the rulers; greater than that! This is because more 
structures come into play at lower length scales; or in other words, the coastline is 
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Figure 1. Fractal structure of Black sea coastline at I§ik (Istanbul), represented actually by a time 
series of stock price. This shows the remarkable self-similarity (examplified by the blow up of the 
segment in the box) involved in many such natural processes. 




Figure 2. Sierpinski gasket (at generation number n = 5), having fractal dimension log 3/ log 2 
(when n — » oo). 

not really a 'line' but rather a fractal having dimension greater than unity. Looking 
at smaller and smaller length-scales, one can find self-similar structures (see Fig. 
!)• 

One can easily construct such regular fractals as carpets or gaskets (see Fig. 
2). In Fig 2, a basic unit of equilateral triangle with an inner triangle obtained 
by joining the mid-points of each side and keeping the inner space void, one 
constructs a gasket. At each step, as the length of each side changes by a factor 
L - 2, the mass of the fractal changes by a factor M{- L d f) = 3, giving therefore 
the fractal dimension of the object to be df - log 3/ log 2. This object in Fig. 2 
of course represents a non-random fractal. However, a random fractal (as in Fig. 
1) can be easily constructed if the void is not always at the center but at any of 
the 4 triangles at random for each generation; the (mass) dimension df remain 
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Figure 3. Construction process of Cantor set is shown (upto n = 3). The set becomes a fractal 
when n — > oo. 

the same. A similar fractal of dimension log 2/ log 3 can be constructed as the 
Cantor set shown in Fig. 3. Starting from a set of all real numbers from to 1 one 
removes the subset from the middle third, so that for each / = 3, the mass (size) 
of the set m - 2 and hence the dimension as n —> oo. This void subset can again 
be randomly chosen, giving a random Cantor set, with same df. 

3. Fractal overlap model of earthquake 

Overlapping fractals form a whole class of models to simulate earthquake dy- 
namics. These models are motivated by the observation that a fault surface, like 
a fractured surface (Chakrabarti and Benguigui, 1997), is a fractal object (Okubo 
andAki, 1987; Scholz and Mandelbrot, 1989; Sahimi, 1993). Consequently a fault 
may be viewed as a pair of overlapping fractals. Fractional Brownian profiles 
have been commonly used as models of fault surfaces (Brown and Scholz, 1985; 
Sahimi, 1993). In that case the dynamics of a fault is represented by one Brownian 
profile drifting on another and each intersection of the two profiles corresponds 
to an earthquake (De Rubeis et al, 1996). However the simplest possible model 
of a fault - from the fractal point of view - was proposed by Chakrabarti and 
Stinchcombe (Chakrabarti and Stinchcombe, 1999). This model is a schematic 
representation of a fault by a pair of dynamically overlapping Cantor sets. It is 
not realistic but, as a system of overlapping fractals, it has the essential feature. 
Since the Cantor set is a fractal with a simple construction procedure, it allows 
us to study in detail the statistics of the overlap of one fractal object on another. 
The two fractal overlap magnitude changes in time as one fractal moves over the 
other. The overlap (magnitude) time series can therefore be studied as a model 
time series of earthquake avalanche dynamics (Carlson et al, 1994). 

The statistics of overlaps between two fractals is not studied much yet, though 
their knowledge is often required in various physical contexts. It has been es- 
tablished recently that since the fractured surfaces have got well-characterized 
self-amne properties, the distribution of the elastic energies released during the 
slips between two fractal surfaces (earthquake events) may follow the overlap dis- 
tribution of two self-similar fractal surfaces (Chakrabarti and Stinchcombe, 1999; 
Pradhan et al, 2003). Chakrabarti and Stinchcombe (Chakrabarti and Stinchcombe, 
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Figure 4. (a) Schematic representations of a portion of the rough surfaces of the earth's crust and 
the supporting (moving) tectonic plate, (b) The one dimensional projection of the surfaces form 
Cantor sets of varying contacts or overlaps (s) as one surface slides over the other. 




time 

Figure 5. For two Cantor sets (as shown in Fig. 4(b); n = 2), one moving uniformly over the other 
(periodic boundary conditions), the total measure of the shaded region contribute to the overlap s; 
the time variation of which is shown (for n = 4). 

1999) had shown analytically by renormalization group calculations that for reg- 
ular fractal overlap (Cantor sets and carpets) the contact area distribution p(s) 
follows a simple power law decay: 

PCs) ~ s^; y=l. (3) 

Study of the time (t) variation of contact area (overlap) s(t) between two well- 
characterized fractals having the same fractal dimension as one fractal moves over 
the other with constant velocity, has revealed some features which can be utilized 
to predict the 'large events' (Pradhan et al, 2004). 
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Bhattacharyya (Bhattacharyya, 2005) has recently studied this overlap distri- 
bution for two Cantor sets with periodic boundary conditions and each having 
dimension log 2/ log 3. It was shown, using exact counting, that if s = 2 n ~ k (n is 
the generation number) then the probability p(s) to get an overlap s is given by a 
binomial distribution (Bhattacharyya, 2005) 

^) = („-ik)(|P(|) 4 ~exp(-r^);^0, (4) 

where r 2 = [| (in - kj\ . Expressing therefore r by log s near the maxima of p(s), 
one can again rewrite (4) as 

/ {logs) 2 ] 

p(s) ~ exp ; n — > oo. (5) 

I n ) 

Noting that p(s)d(\og s) ~ p(s)ds, we find p(s) ~ s~?, y - 1, as in (3) as 
the binomial or Gaussian part becomes a very weak function of s as n — > oo 
(Bhattacharyya et al, 2005). It may be noted that this exponent value y - 1 is 
independent of the dimension of the Cantor sets considered (here log 2/ log 3) or 
for that matter, independent of the fractals employed. It also denotes the general 
validity of (3) even for disordered fractals, as observed numerically (Pradhan et 
al, 2003; Pradhan et al, 2004). 

Identifying the contact area or overlap s between the self-similar (fractal) crust 
and tectonic plate surfaces as the stored elastic energy E released during the slip, 
the distribution (3), of which a derivation is partly indicated here, reduces to the 
Gutenberg-Richter law (2) observed. 

4. Summary 

We introduce here the two-fractal overlap model of earthquake where the average 
distribution of the overlaps between the surfaces, as one fractal (here, Cantor 
set) moves over the other gives the Gutenberg-Richter like distribution (3), with 
y=l exactly in the model. We note that this is an exactly solvable model of 
earthquake dynamics and the result for the distributions compare favorably with 
the observations. 
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